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Semileptonic decays in the limit of a heavy daughter quark 
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The rate of the semileptonic decay b civ is calculated with O (a^) accuracy, as an expansion 
around the limit of equal masses of the b and c quarks. Recent results obtained around the limit of 
the c-quark much lighter than b are confirmed. Details of the new expansion method are described. 



I. INTRODUCTION 



Very recently, next-to-next-to-leading order (NNLO) QCD corrections to the semileptonic h ^ c decay were calcu- 
lated with full account of the charm quark mass [11,12] (see also Ref. Q). The former paper employed a numerical 
method while in the latter the decay rate was expanded in the ratio of quark masses, p = mc/mb. These two ap- 
proaches are complementary, with the numerical one having better accuracy for the larger daughter quark mass and 
04 , the analytical expansion being obviously better for a lighter one. The physically most interesting is the region of the 
actual quark mass ratio, p ^ 0.25.. 0.3. Both methods are applicable in this domain and agree very well with each 
other. The resulting prediction for the 6-quark decay rate will improve the accuracy of the quark mixing parameter 
Kb. 

^SJ \ In the present paper we provide an additional check of that QCD correction. We construct an analytical expansion 
like in [2], but around the opposite limit: instead of starting with p ~ Q (massless charm) we expand around p = 1 
(equally heavy h and c quarks). We find that this leads to a faster convergent series whose sum smoothly matches 
: that found in [2]. As a result we now have a set of analytical expressions valid in the whole range of possible quark 
I masses. In addition, the method of asymptotic expansions is applied to a new kinematic configuration. 
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II. CALCULATIONAL METHOD FOR THE DECAY RATE 

^ , Using the optical theorem, we calculate the decay width from the imaginary parts of 6-quark self-energy diagrams 
' up to four loops, such as shown in Fig. [TJ These diagrams contain two masses, rrib and nic^ and it is not known how to 
\ compute them analytically. We thus treat the mass difference Wf, — rric = '7ib(l — p) as a small quantity and construct 
. an expansion around the limit of equal masses. The expansion parameter is (5 = 1 — p. 

' This expansion is peculiar in the sense that the decay is not possible at the limiting point, ^ = 0. This leads 
to a strong suppression of our result as p tends to one and, as will be seen below, ensures good convergence of the 
expansion. Furthermore, there are no contributions from the region where all loop momenta are of order mf,. This 
OO ', makes the calculation significantly simpler than the complementary expansion around p = 0] ■ 

' A somewhat similar configuration was considered in [4] , where the decay h uiv was evaluated near the limit of 
the maximum invariant mass of the leptons. The difference in the present case is that it is the daughter quark that 
is massive and almost saturates the phase space. Since that massive quark radiates, the calculation is more involved. 
We explain it in some detail below. 

In the first step we integrate over the loop momentum in the massless neutrino-charged lepton loop, replacing it 
with a fractional power of the momentum fiowing through it, l/fc^*^ where Z? = 4 — 2e is the number of dimensions in 
dimensional regularization. For example. Fig. [21(a) shows the lowest-order diagram; the lepton-loop momentum is /. 

The remaining loop momenta can have one of two characteristic scales, hard rrib or soft mb — rric = 5- nib. Depending 
on their configuration, the propagators can be expanded in some small parameter, leading to a factorized product of 
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Figure 1: Sample diagrams contributing to the decay width. Thick and thin lines denote b and c quarks, respectively. Curly 
lines denote gluons and dashed lines denote leptons. All quark flavors have to be considered in the closed loop. 
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Figure 2: In (a) the tree level diagram for the decay b clu is depicted, while (b) shows the general topology of the first 
order correction. Thin and thick lines indicate c and h quarks, respectively. Dashed lines denote leptons and curly lines denote 
gluons. 
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|fc| ~ ruh [1] ^ fc^ + 2po • fc 
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|fc| ~ Srrn, [1] -^2po-k + 2&pI 



Figure 3: The asymptotic expansion of the diagram in Fig. [2] (a) . The dotted double line corresponds to the epsilon propagator 
and the double solid line corresponds to the eikonal propagator 2po • k + 25pq 
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one or more single one-scale integrals [1, H, 0]- This procedure is illustrated in Fig. [3] with the lowest-order example 
without gluons. 

In this example. Region 1 contributes only to the real part, hence need not be considered in this calculation of 
the decay rate. More generally, the hard regions (when all momenta are hard, ~ mf,) will not contribute even at 
higher orders in as . This removes what would otherwise be the most difficult part of the calculation. In the present 
expansion around 6 = there are fewer regions that need to be considered than in the expansion around p = 0. At 
most four regions contribute to a given diagram, compared to eleven in the complementary expansion of Refs. 0,0]. 
All other regions are either purely real or scaleless and give no contribution. 

The second region, thus, contains full information about the tree level result. The resulting integral has already 
been considered in the literature Q. Here we treat its generalization since it will be needed in the higher order 
corrections, 

/d^k 1 
{2tt)D (fc2) Ai (2p„ . fc)A. (2po • k + A)^3 

^ ( 2^A,-g/2 A g-2A,-A.-A, ' 2Ai - A2)r(2Ai + + A3 - D)T{D/2 - Ai) 

^ (4^)^/2r(Ai)r(A3)r(£)-2Ai) ' ^ ' 

where A = 2(5pQ — iO, Pq = ^"mf^, and Ai, A2, and A3 are arbitrary complex numbers. In general, there would also be 
scalar products in the numerator, but it is well known how to deal with these 8] and bring the integral into the form 
of Eq. (IT|). This integral is one of only five general integrals needed in NLO and NNLO calculations. The other ones 
are on-shell propagators up to two loops and one- loop massless propagators, all of which are well known 

For every topology, the most complicated integrals were encountered in the regions where all loop momenta are 
soft. Fortunately, these three- loop integrals could easily be written as a combination of nested integrals of the 
form /(Ai,A2,A3). To illustrate this, let us consider the general two- loop integral, corresponding to Fig. [5] (&) after 
integration of the lepton loop. If both loop momenta are soft, the integral is given by 
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where the are integer numbers and qq is always positive. The ^2 integral can be carried out using Eq. ([T]) with A 
replaced by 2pQ • A:i + A, performing tensor reduction in the process. The resulting ki integral is again of the type of 
Eq. (P). 

For the NNLO calculation, it turned out to be useful to apply partial fraction decomposition in some cases. For 
example, in Eq. ([3]) we could also use the identity 

1 1/1 1 

(2po • k2) [2po • (fci + fca) + A] ^ (2po • fci + A) ^^2^0-^2 " 2po • {ki + fe) + A 

to reduce the number of terms in the denominator. Note that the k2 integral becomes scaleless for 05 < 0. While it 
is obviously not necessary to apply Eq. (|4]) in the case of the integral in Eq. ([3]) , it was necessary to apply analogous 
identities in order to write some of the NNLO integrals as nested integrals of the type of Eq. ([1]) . 

New types of integrals appear only in the diagrams with three-gluon interactions (cf. Fig[T] (c)). due to the third 
gluon propagator. However, in these cases it was possible to apply the so-called Laporta algorithm [l^, [ll| to dispose 
of one of the three gluon propagators. The remaining integrals were again a nested set of /-type integrals. For this 
reduction we used the C++ program rows [12] and the Mathematica package FIRE [3]. 

Our calculation was performed with two independent setups. One approach is based on the code developed for the 
calculation of Ref. [4] . The other uses QGRAF [lj| to generate the diagrams and q2e and exp , to process them 
further (no expansion is done in this step) . The final calculations are in both cases done with custom code written in 
FORM [l||. 




III. RESULTS 



The result for the total width can be cast into the form 

Gj,\Vcb\^ml \ as{mh) 
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where Gf is the Fermi constant and the ellipsis denotes higher order contributions, is defined with five active 
flavors. In QCD we have Cp = 4/3, Ca = 3, and Tp = 1/2. ni — Z denotes the number of light quark flavors, which 
are taken to be massless in our calculation. X^ and Xf, denote the contribution from self-energy diagrams with closed 
c- and 6-quark loops, respectively (cf. Fig[l](a)). Thus, Xc also contains contributions from real c-quark pairs. The 
quark masses are renormalized in the on-shell scheme. 

The tree level and one-loop contributions can be inferred from the closed- form result of Ref. p^ . Expanded in 5 
they read 
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where the ellipses denote higher order terms. Note that the expansion starts at the fifth power of 5. Thus, the total 
width tends to zero very fast as 5 (p) tends to zero (one). Logarithms of S always appear as In 25, since they stem 
solely from A in the integral of Eq. ([T]). 

The first three terms in the expansion of the individual contributions read 
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Figure 4: X2 as a function of p. The solid line denotes the result of Ref. 2|. The dashed lines denote our expansion through 
order (5^ 5^ S^°, and The latter two are almost indistinguishable. (On the vertical axis, zero does not coincide with the 
upper edge of the plot.) 
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We have calculated the fermionic contributions Xi, Xc, and Xi, through terms of order S^^ , while we computed terms 
of order 6^^ and 6^^ for the abelian and non-abelian contributions. Higher order terms are not shown for brevity, but 
are available among the source files of this paper in arXiv. 

To illustrate the convergence behavior of our result, Fig. 2] shows the full NNLO contribution, X2, as a function 
of p. It shows the expansion truncated at different orders in S compared to the result of Ref. [2] . The latter is only 
given up to p = 0.255, which is were the results are closest. The convergence behavior of the expansion around 
p = was studied in Ref. Due to the suppression of our expansion at small values of 5, the different curves are 
indistinguishable for p > 0.4. However, the convergence behavior is very good even close to p = 0. This is in contrast 
to the expansion of Ref. [2], which tends to ±00 as p tends to one. 

Fig. [5] compares the different color structures with the expansion of Refs. 0, Q. As the transition point between 
the two results, we chose the point were they are closest. The two expansions match very well for p between 0.2 and 
0.4. Thus, a combination of the two results enables us to describe the decay over the whole range of kinematically 
allowed values of the daughter quark mass. It was noted in Ref. [3] that the contribution from closed 6-quark loops 
shows an extremum around p — 0.2 (cf. the last panel in Fig. [5]). Using our expansion through <5^^, we were able to 
verify this behavior. 



IV. CONNECTION WITH THE ZERO-RECOIL FORM FACTOR 



In this Section, we provide an independent derivation of the first two terms in the S expansion of the b c decay 
width. They are independent of the real gluon radiation. The real radiation is suppressed by the square of the velocity 
of the daughter quark and influences only the third order terms, of relative O (<5^). The first two terms, O [6^'^), 
are determined by the form factors rjAy describing the M^-boson coupling to quarks. Those form factors arise from 
virtual corrections and replace 1 — 75 in that coupling by rjy — rjAlz- The decay rate expanded in S is, in the lowest 
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Figure 5: NNLO contributions of the different color structures of the total width as functions of p. Thick and thin lines denote 
the results of Eqs. (f9))- (|13|) and Ref. respectively. Note that our expansion tends to zero very fast as p tends to one. (On 
the vertical axis, zero does not coincide with the edge of the plot.) 



two orders, fully described by these form factors, 

, ^ Gl\V.\^mlir^,^,,l) {^)] . (14) 

Without strong interactions, yyy — riA — ^ and we reproduce the first two terms of Eq. ([7]). 

Both form factors are functions of , where q is the four- momentum released in the decay. Thus, to be precise, 
we should have used certain average values in Eq. p4)) . However, when the quark masses are close to each other, the 
variation of is of the second order in 5 and can be neglected in our approximation. 

Even at a fixed q^, the form factors depend on the difference of the quark masses. However, in our approximation 
it is sufficient to know them in the limit of equal quark masses: because of the symmetry to;, ^ TOc, the linear term 
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in S vanishes, and the dependence on the quark mass difference starts only with the quadratic term. In this hmit rjy 
equals one to all orders, while tja is modified by the strong interactions at O (as) and higher orders. Those corrections 
were calculated in Ref . fl^ w ith two-loop accuracy, and in Ref. [13] at three loops. (Full q'^ dependence at two loops 
can be found in Refs. [UUll [H, (H.) 

In order to compare our results with those of Ref. [l^ , we have to change the renormalization scale of as ■ While we 
used as{mb) in Eq. Ref. [l^ uses as{y^mcmt,). Note that a mistake in the running of in Ref. [l^ was pointed 
out in Ref. psj . In the running from nib to ^/rn^rrih, five instead of four active flavors were used. To correct for this 
mistake, we run as in the result of Ref. 19] from y/rn^mh to mi,, using five flavors. At the scale mf,, we decouple the 
b quark and run back to ^JrrQrrvb, using now four active flavors. Thus, the correct result is obtained by adding 
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to 77^ y in Ref. [l^ . In our comparison the correction term contributes to the term of relative order O (S) . 

For completeness we provide all terms of Ref. [l9l | which are needed for the comparison with our result. The QCD 
corrections to the axial form factor are defined as 
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where is defined with four active flavors, rj^ combines the contributions from diagrams with closed c- and 6-quark 
loops. The individual color structures in the limit S ^ are given by 
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The term of order 6 in rj^ is due to the correction term in Eq. (jlSp . This linear term arises because the m(, 
symmetry is broken by the charge renormalization, since the &-quark does not contribute to the running between rrib 
and y/m^jrib. 

To compare the two results, we decouple the b quark in our result and run as from mf, to y/m^Jfib with four active 
flavors. This changes X2 by 
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We find perfect agreement for the first two terms of our expansion. Comparing the widths calculated with Eqs. ([5]) 
and (|14p . we expect and indeed confirm that 
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In the tree-level term Xq, the factor 3/4 eliminates the contribution of the vector coupling. 

Individual parts in Eqs. ((9l ll3p are also reproduced; r]^ combines effects of both heavy quarks and corresponds to 
the sum of Xc and Xb . 



V. SUMMARY 



To summarize, we have calculated the semileptonic b c decay as an expansion around the limit of equal quark 
masses. Our result is a fast convergent series, which smoothly matches the expansion in the opposite limit. Our result 
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confirms the calculations of Refs. [l!, [2]. Together with the result of Ref. 0, we now have analytical results valid 
over the whole range of kinematically allowed daughter quark masses. An additional check of a part of our result is 
afforded by comparing with the result of Ref. [l^ . 

Furthermore, we have explained the application of the method of asymptotic expansion to a new kinematic limit. 
This limit leads to significant calculational simplifications and results in a fast convergent series, which is applicable 
over most of the allowed region of the daughter quark mass. Even in the massless limit the error is only about 2.5% 
for the total width, as demonstrated in Fig. [J] Thus, this new expansion provides a convenient tool for future studies 
of various aspects of decays not only of quarks, but also leptons such as the muon. 
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